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The formation of a core collapse supernovae (SNe) results in a fast (but non- or mildly-relativistic) 
shock wave expanding outwards into the surrounding medium. The medium itself is likely modi¬ 
fied due to the stellar mass-loss from the massive star progenitor, which may be Wolf-Rayet stars 
(for Type Ib/c SNe), red supergiant stars (for type IIP and perhaps Ilb and IIL SNe), or some other 
stellar type. The wind mass-loss parameters determine the density structure of the surrounding 
medium. Combined with the velocity of the SN shock wave, this regulates the shock acceleration 
process. In this article we discuss the essential parameters that control the particle acceleration 
and gamma-ray emission in SNe, with particular reference to the Type lib SN 1993J. The shock 
wave expanding into the high density medium leads to fast particle acceleration, giving rise to 
rapidly-growing plasma instabilities driven by the acceleration process itself. The instabilities 
grow over intraday timescales. This growth, combined with the interplay of non-linear processes, 
results in the amplification of the magnetic field at the shock front, which can adequately account 
for the magnetic field strengths deduced from radio monitoring of the source. The maximum 
particle energy can reach, and perhaps exceed, 1 PeV, depending on the dominant instability. The 
gamma-ray signal is found to be heavily absorbed by pair production process during the first 
week after the outburst. We derive the time dependent particle spectra and associated hadronic 
signatures of secondary particles (gamma-ray, leptons and neutrinos) arising from proton proton 
interactions. We find that the Cherenkov Telescope Array (CTA) should be able to detect objects 
like SN 1993J above 1 TeV. We predict a low neutrino flux above 10 TeV, implying a detectability 
horizon with current or planned neutrino telescopes of 1 Mpc. 
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1. Introduction 


Supernova Remnants (SNRs) are widely accepted as the likely source of Galactic cosmic rays 
(CR) up to energies of 3 PeV (the CR “knee” §)■ Per the Hillas confinement criterion [Q] that 
the Larmor radius of the particle match the source size, only certain Galactic sources should be 
able to produce such energetic particles. These include extended sources with standard interstellar 
medium (ISM) magnetic field (MF) values such as massive star clusters and their superbubbles 
till; or compact sources with more intense magnetic fields, such as young SNRs. In this 
work we explore particle acceleration and gamma-ray radiation in SNe arising from massive star 
progenitors, with particular focus on the well-monitored SN 1993J. We consider the acceleration 
of particles starting from the days following SN outburst when the radio luminosity is close to its 
maximum. We review the main results of radio observations and modeling in §0, examine models 


of magnetic field amplification (MFA) in SN 1993J in §2.1, and the implied cosmic ray energies 


in §^. Modeling of different particle distribution species in radio SNe is outlined in § |2.3| . Our 
calculation of the gamma-ray radiation is shown in §|3|, and the neutrino signal in §3.3. Some 
perspectives for other types of SNe, along with conclusions, are discussed in §^ 


2. A case study: SN 1993J 

Our model for SN 1993J is based on work by [^(hereafter T09) that discusses particle accel¬ 
eration in SN 1993J. This Type Ilb SN resulted from the explosion of a massive K-supergiant star 
with an initial mass ~ 15 Mq 0, having a B-star binary companion 0. The progenitor star had 
a mass loss rate ~ 3.8 x lO^^M© yr^^ and a wind velocity u„ ~ 10 km s^^ (T09). Assuming a 
constant mass loss rate, the circumstellar medium (GSM) density scales as tidrc = — 


4;rr-WH’M//(l+4X) 

where A = 0.1 is the Helium fraction and mu the hydrogen atom mass. The effective density 
downstream of the forward shock, one day after outburst, is ~ 4 x 10^ cm^^ with a shock 
compression ratio of 4. The shock radius after 1 day has been estimated as r ~ 3.5 x lO'^ cm (T09) 
from the radio expansion. The shock propagates into a fully ionized medium 0. 

The time dependence of the average MF from the synchrotron emitting shell is (T09): 


{B) ~ [2.4 ± 1 G] X 


- 1 . 16 ± 0.20 


( 2 . 1 ) 


100 days^ 

giving a MF of order 500 G after 1 day. This field is consistent with observations [ |T^ ] . It is an 
averaged value over the synchrotron shell and does not necessarily represent the MF produced at 
the forward shock. We identify this field as that in the post shock gas of the forward shock front. 

2.1 Magnetic field amplification (MFA) 

Eq. ( ^!1| ) shows that the MF substantially exceeds typical stellar wind fields at these distances. 
For comparison the equipartition MF in the progenitor wind is 0: 


Beq = 


[2.5 mG] X X u^ jg x r^g' , 


( 2 . 2 ) 


about one thousand times less than the value in Eq.^. T09 find fhat the sub-shock compression 
ratio is close to 4. The compression of the ME given by Eq. (p.2[) cannot therefore explain the value 
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given by Eq. (2.1) (see also [pT]]). We conclude and assume that a strong MF amplification process 
is at work at the forward shock front, driven by the diffusive shock acceleration (DSA) of hadrons. 

2.1.1 Streaming driven instabilities 

Streaming of cosmic rays ahead of the shock front produces magnetic fluctuations. The stream¬ 
ing modes can be in resonance (R) with the energetic particles, i.e. in the high-energy limit they 
have a wave-number such that k ~ rT ^; or they can be non-resonant (NR) with a much larger wave- 


number [12, 13]. The NR modes grow the fastest 114]. For 1993J the NR modes grow and produce 
magnetic fluctuations over intra-day timescales; we obtain a growth time: 


TlNR-st = [0.16 day] x 


0/15 




CR 


/0.05)u3 


sh.93J 




X EpeV tdayl 


(2.3) 


The CR distribution is assumed to scale as p^^ over more than 6 orders of magnitude producing 0 = 
^t\-{pmax/Pinj) = 15, and ^cR = 0.05, so that 5% of the fluid kinetic energy is imparted to energetic 
particles. p„un: (Pinj) is the maximum (injected) particle momentum. The growth timescale has to 
be shorter than the advection timescale towards the shock Zadv = latter is calculated for 


a diffusion coefficient K = rj Kg > taken in the background MF (Eq ^ ) and Kg = cr^f 3. We 
find: Tadv = [0.24 day] x rj Epev x The condition X^Rst < '^adv is necessary but not sufficient 
for the instability to develop. The magnetic fluctuations produced by the instability uncovered by 


[ ]15| ] produce small scale perturbations. The wave number corresponding to the maximum growth 
rate is [ 13] komoxfimax — 4 X 10^ for SN 1993J. Acceleration and confinement of energetic particles 


up to a few PeV requires magnetic fluctuations to be generated at resonant scales kri^ ~ 1, as the 
NR growth rate scales as the MF at the scale of interest grows over a timescales about 2000 
times larger than the one obtained in Eq. If the background MF is purely toroidal rj < 1, and 
the instability may not have time to develop. 

The R instability can build up over the NR one [jf^, with the ratio of the magnetic ener¬ 
gies reaching s/^^^cjv^h 0.95 in our case. The growth timescale for the R instability is however 
longer than the rate given by Eq. 2.3. We have [ |T3| ] XRst — sjTto!% / r^max with a ~ 3 x lO^^cm^/s^ 
in the conditions that prevail for SN 1993J. This leads to a growth rate Xr-si — 16 x^Rst at 1 
PeV. The NR instability can be driven to non-linear stages and produce large wave numbers. A 
typical timescale of non-linear saturation of the MF is about 5 x X^R^st Recently [ jT^ ] pro¬ 
posed a ponderomotive instability that builds up on the magnetic fluctuations by the NR stream¬ 
ing instability. We can evaluate the growth time-scale of such long-wavelength modes: Xrw = 

^ X Epev tdays- The parameter A = > 1 is the level of 




\/“sii,9: 


jAio 


magnetic energy produced by the small scale instability with respect to the background CSM and is 
normalized to 10. Fong wavelengths are produced on timescales shorter than Xadv for kri^max — 1- 


2.1.2 Turbulence driven instabilities 

Stellar winds of massive stars are subject to strong fluid instabilities that can lead to turbulent 
motions [j^. Turbulent density and magnetic fluctuations at SN shocks can result in MF amplifica¬ 


tion [j^, 20]. In the latter work the presence of CRs is not necessary to produce the magnetic field 
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amplification. The MF growth time is controlled by the coherence length of the turbulent spectrum 
L and the fluid velocity Ush- For SN1993J, L/ush [0.3 year] Lo.oipc/ush,93J hence the magnetic 
field has to grow over a fraction of 10^^ of this timescale. 


2.2 Maximum cosmic ray energies 


Immediately after SN outburst, the maximum energy is likely limited by the SNR age. Bal¬ 
ancing the age with the acceleration time gives Tacc = g{r)i^u/u^sh’ "'hh r the shock compression 
ratio, and g{r) = 3rl{r — 1) x (K'd/K:„r+ 1). We use Kd = JCu/v/TT, corresponding to a tangled 
MF whose tangential component is compressed by a factor 4 111]- The maximum energy is 
Emax,age,PeV ^ X (1 But rapidly (see the streaming instability amplifies the 

MF, and the non-linear process produces a MF at saturation. The typical saturation value [|T^] is: 


= [16 Gauss] x ^ x t^^ . (2.4) 

This value is within a factor of 2 of the MF derived in the upstream medium from Eq.^T] using a 
compression ratio r = 4. If only the NR instability is at work building the MF, the maximum parti¬ 
cle energy is then fixed by a condition over the CR areal charge [ |I^ that produces / = 6.8. 

The latter value corresponds to the amplification of the equipartition MF to the value deduced from 
radio observations. In that case: Emax,NR,PeV ~ 1 x In the case of long wavelengths fluc¬ 

tuations produced by the ponderomotive instability, the maximum energy is fixed by geometrical 
losses. The maximum energies are obtained with a diffusion coefficient expressed in the amplified 
field and compared to T\escEshitsh- In order to derive a time dependence of the maximum energy the 
time dependence of the amplified field has to be specified, for which we rely on the estimate given 
in Eq. This gives: Emax,LW,PeV ~ 55 (^) x ■ Apart from the time dependence of Bsat, 
this value is optimistic since the highest energy particles may feel a ME lower than B^at [^. The 
maximum cosmic ray (proton) energy at any given time is the minimum of all the above limits. In 
all cases, PeV energies are possible days after the outburst. 


2.3 Cosmic-ray spectral evolution 

We follow the prescription of T09 regarding the time evolution of the CR energy content. 
The proton particle spectrum is assumed to follow a power-law with a spectral index s = 2 and an 
exponential cutoff at Emax{t), the maximum particle energy as discussed in §^. In the case of 
secondary electrons and positrons produced in the proton-proton interactions, we solve a one-zone 
energy equation to calculate the time evolution of their energy distribution N{E,t): dtN{E,t) -Y 
dE{L{E)N{E,t)) = Q{E,t), where L includes the synchrotron losses for secondary electrons and 
positrons in the post-shock region. 


3. Gamma-Ray production 

We derive multi-wavelength time dependent spectra, in the context of MEA, driven by the 
instabilities discussed in section § |2.1[ We consider only proton-proton interactions, as inverse 
Compton or bremsstrahlung radiation have been found negligible in the GeV-TeV range explored 
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here. Inverse Compton process is highly disfavored due to the strong magnetic field at the forward 
shock. Gamma-rays can be absorbed by different soft photon fields to produce electron-positron 
pairs. The main photon source is the SN photosphere, described in the case of SN 1993J by |^]. 


3.1 Gamma-gamma absorption 

We have performed a full calculation of the gamma-gamma opacity Tyj including geometrical 
effects due to the anisotropic interaction (Renaud et al. 2015, in prep). The final gamma-ray flux is 
the unabsorbed flux TV,™ times an attenuation factor exp(—Tyy(£'y)). Gamma-gamma absorption is 
strong just after the outburst as the interaction in nearly isotropic, but thereafter drops as the ratio 
of the forward shock radius to the photosphere radius reaches ~ 3, which happens after ~ 5 days. 


3.2 Cherenkov Telescope Array detectability 

The time dependent gamma-ray spectra in the very-high energy (100 GeV < E < 100 TeV) 
gamma-ray domain are displayed in Figure 1. A source like SN 1993J would be easily detected by 
CTA above 1 TeV in 20 h of observing time. The best time window to detect a gamma-ray signal 
is between a week and a month after the outburst. Prior to a week the source is optically thick 
to gamma-rays, but the gamma-gamma opacity decreases rapidly due to anisotropic effects. After 
a month the gamma-ray signal becomes too faint due to the decrease in density of the medium. 




Figure 1: [left] Time dependent spectra at 4 different times after outburst. Dotted lines - expected CTA 
sensitivity in 50h. At t= 1.18 days, the source is optically thick to gamma rays, and the flux is below the given 
range. [Right] Time dependent neuttinos flux above an energy E expected by a KM3NeT-like instrument 
(continuous lines) at 4 different times after outburst. Dotted lines - atmospheric neutrino backgrounds. 


3.3 Neutrino signal 

Neutrinos are by-products of pion production. The expected flux of neutrinos detected by an 
instrument equivalent to KM3NeT is displayed in Figure 1. We find that at best about 0.1 neutrinos 
could be expected from a SN 1993J type event above 10 TeV by summing the spectra between 1 
and 30.28 days after the outburst, with 0.03 background neutrinos. A signal of one neutrino from 
such a source requires a source within ~ 1 Mpc, or a gamma-ray signal ten times stronger. 
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4. Discussion & Conclusions 


SNe produce multi-PeV particles through the combination of: fast shocks (v ~ 0.05c), high 
density CSM produced by stellar winds, and low wind magnetizations. High degree of CSM ion¬ 
ization eases the particle acceleration process. Assuming that the background MF has a turbu¬ 
lent component, different instabilities driven by the acceleration process can grow over intra-day 
timescales. This model is consistent with the strength of the MF seen in SN 1993J. 

Parameters affecting the early gamma-ray emission from SNe include the ratio of the mass 
loss rate to the wind velocity {M/v„), which fixes the CSM medium density and affects the CR 
driven instability growth rate. The shock velocity controls the growth rate of the instabilities and 


the acceleration timescale. The degree of ionization [ 24 ] is important for the particle acceleration 
efficiency and may also produce elemenf dependenf CR specfra if fhe ionizafion is partial. The 
background MF confrols parfly fhe local magnefizafion and fhe shock obliquity. The SN luminosify 
confrols fhe gamma-gamma absorpfion process. 

Only abouf 5-6% of fhe local core-collapse SNe have been classified as Type Ilb such as SN 
1993J [S]. There appear fo be fwo classes of lib SNe, wifh compacf and exfended progenifors 


[]26|]. If is likely fhaf only fhe ones wifh extended radii and higher mass-loss rates such as 1993J are 
fhe more likely candidates for defecting y-ray emission, fhus furfher reducing fhe sample. Type Iln 
SNe are probably more promising fargefs for gamma-ray telescopes in ferms of high ambienf den¬ 
sity wifhouf significanfly reduced velocities. In fhe case of fhe Iln SN 1996cr, if has been deduced 
fhaf fhe shock was inferacfing wifh a shell of density ~ 10^ cm^^ |27] a few years afler explosion. 
This provides a high densify largel for producing y-rays via pion production (allhough perhaps nol 


as dense as suggested by some aulhors [ 128| ] for olher SNe). Such Type Iln SNe would probably 
also be promising fargefs for defecting neufrino emission from secondaries. Unforlunafely fhese 
sources are even less frequenf fhan fhe lib SNe. SN IIP comprise fhe largesl class of core-collapse 
SNe, making up around half fhe folal. They arise from RSG slars, which have wind mass-loss 
rates ranging from 10^^ fo 10^“^ Mq yr^' [|2^. However, observafionally mosf IIPs appear fo arise 
from fhe lower end of RSG sfars [ ^5| ] and are less luminous in X-rays [M|, suggesfing fhaf only 
fhe RSG slars wifh initial masses below ~ 17-19 Mq, wifh correspondingly lower mass-loss rates 


], explode fo become Type IIP SNe. The rare Ib/Ic SNe harbor fhe faslesl shock waves, bul arise 
from Wolf-Rayel progenifors, which have wind velocifies fwo orders of magnilude greater fhan 
RSGs, and fherefore should have a correspondingly lower wind densify. On fhe ofher hand, Iheir 


X-ray flux is presumed fo be due fo Inverse Compfon or synchroiron emission [30|, suggesfing 
accelerated elections, so if is possible fhaf fhe shocks are capable of also accelerating prolons fo 
high energies. W-R slars are surrounded by low density wind-blown bubbles bordered by a high 
densify shell. If fhe shell is formed soon before fhe explosion, as is fhe case of fhe SN 2006jc [^], 
Ihen if provides a good largel for acceleraled prolons fo collide wifh. Such W-R slars may be good 
candidales for delecling gamma-ray emission in fhe early phases. There may be SNe similar fo SN 


1987A, whose progenifor, a blue-supergianl, had a very low mass-loss rate wind ||32[], bul which 
shows evidence for a dense HII region wifh density of order 200 parlicles cm^^ [ 331, surrounded by 
a dense circumslellar ring wifh density ~ lO'^ particles cm^^ al a dislance of ~ 0.2 pc from fhe SN. 


Finally, fhe class of super-luminous SNe, especially Ihose fhaf are H-rich [34], may be inferacfing 
wifh extiemely dense environmenfs. High densilies close in fo fhe sfar could again provide large! 
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material for proton-proton collisions and detectable 7 -ray emission at an early age. 
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